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The taxonomic status of the central Peruvian Telmatobiinae Is reassessed 
by analyzing the intraspecific variation of 18 morphometric measures among 
the currently recognized taxa (3 genera, 6 species, 9 subspecies), Cluster 
analysis, principal component analysis and discriminant analysis lead to the 
recognition of two genera (Batrachophrynus and Telmatobius) including six 
species (Batrachophrynus brachydactylus, B. macrostomus, Telmatobius bre- 
virostris, T. carriliae, T. jelskii, T. rimac) without segregation in subspecies. 
Two diagnostic external features distinguish Batrachophrynus from Teimato- 
bius species, another two characters are convergent adaptations distin- 
guishing the stream-inhabiting ecotype from the lake-inhabiting one. 


INTRODUCTION 


Leptodactylid frogs of the genus Telmatobius Wiegmann, 1835 inhabit the whole 
range of the Andes from Ecuador in the north to Chile/Argentina in the south (FROST, 
1985; Cet, 1986). However, in the high-Andean regions of central Perú around Lake Junin 
(FseLpsa, 1983), two endemic Telmatobius-like species have been described as members of 
the genus Batrachophrynus Peters, 1873: the large lake-inhabiting Junin frog B. macro- 
stomus and the stout stream-inhabrting B. brachydactylus. PETERS's (1873) distinction was 
based solely on the presence (Telmatobius) or absence (Batrachophrynus) of maxillary and 
prevomerine teeth. Yet, T. brevipalmatus, T. edentatus and T. intermedius lack maxillary 
teeth (LvNcH, 1971). Differences in several osteological characters other than teeth 
between Batrachophrynus and Telmatobius confirmed the generic distinction (LYNCH, 
1978) and indicated an early separation of Batrachophrynus from the Telmatobius stock 
(Cri 1986). LAURENT (1983) assigned B brachydactylus to a monotypic third genus 
Lynchophrys based on morphometric differences from B. macrostomus. At present, this 
view on taxonomy ts maintained, though independent support is absent (DUELLMAN, 1993; 
Frost, 1985; LaviLLA, 1988a). According to LAURENT's opinion, Lynchophrys is more 
recently derived from the Telmatobius stock than is Batrachophrynus. Thus, the 
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phylogenetic relationship among the two species B. brachydactylus and B. macrostomus 
and of both with the genus Telmatobius are controversial. 


The central Peruvian Telmatobius presently include four stream-inhabiting species, T. 
brevirostris Vellard, 1955, T carrillae Morales, 1988, T. jelskii (Peters, 1873) and T. rmac 
Schmidt, 1954, whereas T. juninensis (Shreve, 1938) was an incorrectly classified Phrynopus 
(DUELLMAN, 1993). VELLARD (1955) proposed three subspecies for T. brevirostris, four for 
T. jelskii and two for T. rimac, although the significance of the morphometric differences 
between the proposed subspecies was not tested for, and the existence of intermediate 
individuals was stated. Field studies in the Mantaro Valley (Department of Junín, Perü) 
on T. jelskii raised doubts on the subspecific classification because morphs pertaining to 
different “subspecies” were found within the same population (SiNscH, 1985, 1986, 1990). 


Part of the taxonomic confusion 1n this group is probably due to morphological 
convergence among the species of each ecotype: the stream-inhabiting, semiaquatic frogs 
are stout and usually moderate-sized or small, whereas the lake-inhabiting, fully aquatic 
frogs are large (SINSCH, 1986, 1990). Therefore, phenotypic similarity between different 
species may reflect convergent adaptation to the same environmental factor rather than 
phylogenetic relationship. On the other hand, subtle differences among populations of the 
same species could result from direct responses to such factors as the amount of water flow 
in different streams. To test for environmental influences on morphological traits, we 
complemented our analysis of the two Batrachophrynus species and four stream-inhabiting 
Telmatobius species from central] Perú with the lake-inhabiting Titicaca frog Telmatobius 
culeus (Garman, 1875) from southern Pert. 


It is obvious that the validity of the central Peruvian taxa of Telmatobiinae needs to 
be confirmed by a thorough investigation. The definition of most taxa is exclusively based 
on morphological and morphometric characters (e.g. VELLARD, 1951, 1953, 1955) without 
an appropriate assessment of the intraspecific variation (TRuEB, 1979). Currently, the 
identification of most species 1s only possible by comparison with type specimens 
Consequently, our investigation evaluates VELLARD's type material in MHNSM (Museum 
of Natural History “Javier Prado", Lima, Perú) in comparison with the preserved 
specimens in URP (Museum of Natural History of the Ricardo Palma University, Lima, 
Perü) and with specimens collected in the field in 1992. We use cluster analyses to detect 
intraspecific morphological traits within the data sets of phenotypes described by 
morphometric characters. At the species level we widely follow the procedures used by 
Wiens (1993) in his recent taxonomic revision of the Telmatobiinae from northern Perú. 
Principal component and discriminant analyses are applied to quantify the differentiation 
among the seven nominal species. Finally, the differences among the genera and the 
ecotypes (stream- versus lake-inhabitants) are surveyed in discriminant analyses. The aims 
of our study are to provide an objective basis for the identification of central Peruvian 
Telmatobunae and to distinguish convergent morphological traits from morphometric 
features which characterize species and;or genera. 
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Fig. 1 Batrachophrynus brachydaetylus (A, dorsal view; E, ventral view); B. macrostomus (B, F), 
Telmatobius jelskii (C, Gy, T. rimac (D, H). 
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MATERIAL AND METHODS 


The material examined 1ncluded a total of 280 adult frogs pertaining to the seven 
nominal species Batrachophrynus brachydactylus (figs. LA, 1E), B. macrostomus (figs. 1B, 
IF), Telmatobius brevirostris, T. carrillae, T. culeus, T jelskii (figs. IC, 1G) and T. rimac 
(figs. 1D, 1H). The frogs had been collected at 32 localities (fig. 2). The assignment of 
specimens to subspecies and to localities is listed in Appendix I, the number of individuals 
studied and the sex-specific size are given in Table I. 


Standard morphometric measurements (nearest 0.1 mm) include: (1) snout-vent length 
(SVL); (2) height of body at the pectoral girdle (BH); (3) maximum width of head 
(HWID); (4) eye diameter (EYE); (5) interorbital distance (IOD); (6) eye-nostril distance, 
from anterior margin of eye to posterior edge of naris (ENOSE); (7) distance between the 
eye and the tip of the snout (ESNOUT); (8) humerus length (HUML); (9) radioulnar 
length (RADL); (10) hand length (HNDL); (11) length of the third finger (FG3L); (12) 
femur length (FEML); (13) tibia length (TIBL); (14) foot length, from union with tibia to 
the tip of fourth toe (FOOTL); (15) length of first toe (TOEIL); (16) length of fourth toe 
(TOE4AL); (17) length of callus internus (CIL); (18) maximum length of toe web (WEBL). 


Muluvariate analyses were performed on log, -transformed data (BOOKSTEIN et al., 
1985) and morphometric ratios The empiric measurements were transformed to ratios 
(range: 0-1) by calculating measures relative to SVL (SCHNEIDER et al., 1992, 1993). 
Moreover, two indices were used for further analysis: CIL/TOEIL and FEML/TIBL. 
Table II gives the means (and corresponding standard deviation) of these 19 relative 
measures for each species. 


Due to the low number of individuals assigned to type material, we analysed the 
morphometric similarity between individuals in a cluster analysis using hierarchical 
grouping of the SVL-standardized phenotypes (WARD, 1963). This procedure subsequently 
reduces the number of groups by joining that specimen to another one or to a cluster 
which onginates the lowest error sum of square. The result ıs a dendrogram based on 
phenetic similarity. Ás the proposals of subspecies by VELLARD were exclusively based on 
external morphology, valid subspecific taxa are expected to form homogeneous clusters. 


At the species level, sets of the log,,-transformed data were subjected to principal 
component analysis to explore the morphometric variability independent of taxonomic 
assignment. Principal components (PC) are linear combinations of the measured variables, 
uncorrelated with each other and explaimng the maximum amount of variation. The first 
principal component (PC!) of morphometric data generally describes differences in size, 
but size effects may be present in subsequent principal components (HUMPHRIES et al , 
1981). Techniques such as shearing have been developed to correct PC2 and PC3 for 
possible size effects (BOOKSTEIN et al., 1985), but they are controversial and size effects may 
still persist (ROHLF & BooksTEIN, 1987). Moreover, we applied canonical discriminant 
analysis to distinguish between the taxonomic groups delimited a priori. The resulting 
discriminant functions (CAN = canonical variables) are linear combinations of the 
measured variables that maximize the differences between the groups Discriminant 
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Fig. 2. — Distribution of the central Peruvian Telmatobunae’ Butrachophrynus brachydactylus 


(inverted triangle), B macrostomus (open triangle), Telmatobtus brevirostris (x ), T. carrilae (+), 
T jelskii (circles), F rimac (dots) Localities are approximated from distances by roads, multiple 
localities in close proximity are represented by a single symbol. The main Andean river systems 
are indicated. 
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Tabie |! Morphometric data for Batrachophrynus bracivdacntus, B macrostomus Telmatobius brevirostris T. carriliae T culeus, T jelsku and T 
rimac The first line i$ mean + 1 SD; second Lne i5 range All values are in millimeters see text for abbreviations of variables 


B brachydactylus B. macrostomus T culeus T jeiskat 


T— 
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Table II - Ratios of morphometric data for Barrachophrynus brachydactylus, B. macrostomus, Telmatobius brevirostris, T. carrillae, T. culeus, T 
Jelsku and T, rimac. Data are given as mean + 1 SD See text for abbreviations of variables. 


B. brachydactylus T brevirostris T carnilae T, jelsku 


BH/SVL 0.194 + 0.047 | 0.176 + 0.025 | 0285 + 0023 | 0.221 + 0020 | 0292 + 0032 | 0.252 + 0028 | 0.240 + 0.029 


HWID/SVL 0314 + 0026 | 0.404 + 0.024 | 0.347 + 0.034 | 0.302 + 0.021 | 0.397 + 0.037 | 0.364 + 0.026 | 0.337 + 0.026 g 
EYE/SVL 0.085 + 0.014 | 0060 +0007 | 0103 + 0.004 | 0.085 + 0.012 | 0076 + 0.009 | 0.100 4 0015 | 0106 + 0012 5 
IOD/SVL 0.240 + 0.021 | 0.238 + 0.010 | 0281 +0011 | 0.242 + 0.019 | 0.240 +0014 | 0.270 +0026 | 0.282 + 0.022 m 

ENOSE/SVL | 0.151 + 0.015 | 0 143 + 0.005 | 0.167 + 0.006 | 0.152 + 0.022 | 0.138 + 0.012 | 0.167 c 0014 | 0170 x 0013 v 

ESNOUT/SVL | 0.206 + 0.016 | 0212 + 0.007 | 0236 x 0013 | 0.210 + 0.014 | 0211 £ 0018 | 0.228 +0015 | 0243 + 0.016 m 

HUML/SVL | 0.262 + 0.034 | 0 373 + 0.020 | 0.299 + 0.027 | 0.233 + 0.037 | 0.306 + 0.025 | 0.285 + 0.043 | 0.280 + 0.028 a 

RADL/SVL 0.195 + 0.022 | 0.262 + 0.017 | 0.227 + 0.014 | 0.215 + 0.017 | 0227 £ 0015 | 0.226+ 0023 | 0218 + 0.021 £g 

HNDL/SVL 0.199 + 0.021 | 0.299 + 0.017 | 0.270 + 0.020 | 0217 + 0.027 | 0.257 + 0.017 | 0.255 + 0.023 | 0.260 + 0.022 o 
FG3L/SVL 0.113 +0015 | 0.189 + 0017 | 0.176 + 0.012 | 0128 x 0.017 | 0158 + 0019 | 0.163 + 0.025 | ‘0.166 + 0.024 A 

FEML/SVL 0 440 4 0038 | 0.481 + 0.039 | 0.468 + 0.062 | 0445 + 0.037 | 0.452 + 0.032 | 0.453 + 0.043 | 0.475 + 0.039 > 
TIBL/SVL 0410 + 0035 | 0.441 +0019 | 0.475 + 0.021 | 0.457 + 0.027 | 0467 + 0034 | 0.448 + 0.033 | 0.495 + 0.037 m 

FOOTL/SVL | 0640 1 0047 | 0.729 + 0.039 | 0.741 + 0.051 | 0 713 + 0.042 | 0.751 + 0 062 | 0.695 + 0.070 | 0.763 + 0.044 

TOEIL/SVL | 008840010 | 0.127 £ 0010 | 0.110 x 0.011 | 0.104 + 0.013 | 0108 £ 0013 | 0.089 + 0.011 | 0.117 + 0.0/2 

TOEAL/SVL | 0.418 + 0.030 | 0.458 + 0.033 | 0.489 + 0.047 | 0 469 + 0.029 | 0.506 + 0042 | 0.454 + 0.041 | 0.519 + 0.047 
CIL/SVL 0041 + 0009 | 0035 + 0006 | 0060 + 0.009 | 0.049 + 0.006 | 0047 + 0.007 | 0.047 + 0.011 | 0.051 + 0 008 

WEBL/SVL 0.126 + 0.026 | 0.214 0027 | 0101 + 0.053 | 0 159 + 0.027 | 0.198 + 0 036 | 0.187 + 0.072 | 0.191 + 0.087 

CIL/TOEIL 0.472 + 0.107 | 0.275 + 0.050 | 0.555 + 0.123 | 0 484 + 0.092 | 0.440 + 0.066 | 0.536 + 0.144 | 0.442 + 0 056 

FEML/TIBL 1.072 + 0.074 | 1092 + 0084 | 1 002 + 0.090 | 0,975 + 0.075 | 0.972 + 0.074 | 1.013+ 0.080 | 0.960 + 0 058 

2 


Source . MNHN, Paris 


22 ALYTES 13 (1) 


functions were derived from the log,,-transformed data and from the ratios. The degree 
of separation of taxa was almost identical in both analyses and therefore we present the 
results using the log-transformed data only (analogous to Wiens, 1993). 


All calculations were performed on a PC using the FORTRANT77 program 
CLUSTER and the program package STATGRAPHICS, version 5 5. 


RESULTS 


INTRASPECIFIC VARIATION 


Among the seven nominal species studied in this paper, four are thought to segregate 
in subspecies: Telmatobius brevirostris, T. euleus, T jelskii and T. rimac (VELLARD, 1951, 
1953, 1955). We approached the problem of morphometric variation within a nominal 
species by applying cluster analysis on the phenotypes of all conspecific individuals 
available, including those which VELLARD assigned as type material for the proposed 
subspecies. Phenotypes are described by 19 SVL-standardized morphometric ratios to 
minimize size effects. 


The phenograms of all species (including those without subspecific segregation) 
revealed existence of different intraspecific morphological traits documented by groups of 
specimens which joined to the same cluster with an error sum of squares of less than 0.1 
(figs. 3-6, data on Batrachophrynus and T. carrillae not shown). We treat these groups of 
remarkably similar specimens as “morphs”, without intending a taxonomic implication. 
Generally, intraspecific morphs did not reflect polymorphism related to sex or size. 


Telmatobius brevirostris Vellard, 1955 


VELLARD (1955) recognized three subspecies: Telmatobius b. brevirostris, T. b. parvulus 
aud T. b. punctatus. Unfortunately, we only found in the MHNSM collection two 
specimens of the first two subspecies, and one of the third. The original descriptions were 
based on only three adults of T. b. brevirostris, six of T b parvulus and two of T. b. 
punctatus — a prohibitively small sample size to define reliable subspecific taxa. The 
phenogram (fig. 3) of the five specimens available for analysis shows that: (1) the two type 
specimens of 7. b. brevirostris from Chasqui join different clusters; (2) the morphometri- 
cally most similar specimens are one T b brevirostris individual and the T. b. punctatus 
male from Santa Maria de! Valle; (3) the two individuals of T 5. parvulus from Caina are 
very alike, but form together with one T. b. brevirostris a main cluster with a error sum 
of squares far below the level of morph distinction. All specimens Join to one group at an 
error sum of squares of only 0.146. Neither the association of the five phenotypes to 
groups nor the degree of morphological differentiation between them support a subspecific 
segregation within the nominal species 7. brevirostris. 
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Fig. 3 — Phenogram of morphometric similanty between adult Telmatobius brevirostris which 
VELLARD (1955) assigned to different subspecies Each specimen is identified by the locality of 
collection and an individual number Similarity 1s based on 19 ratios and computed by 
hierarchical grouping in a cluster analyis. T. b. brevirostris: Chasqui, T b. parvulus: Caina, 
T. b. punctatus. Santa Maria del Valle. 


Telmatobius culeus (Garman, 1875) 


VELLARD (1953) recognized six subspecies: Telmatobius c. culeus, T. c dispar 
(redescribed by LaviLLA, 1988b), T. c. escomeli, T. c exsul, T c. fluviatilis and T c 
lacustris. The phenogram (fig. 4) calculated for 42 adults assigned as type material 
(MHNSM ) reveals that: (1) there exist five morphs within the nominal species; (2) these 
morphs do not coincide with any of the subspecies proposed by VELLARD; (3) different 
morphs occur at the same locality. Joining all specimens to one group causes an error sum 
of square of 0.453, that is three times greater than in T brevirostris. Again, from the 
morphometric point of view there is no evidence that the morphological traits within 7. 
culeus agree with the proposed subspecific differentiation. 


Telmatobius jelskii (Peters, 1873) 


VELLARD (1955) recognized four subspecies. Telmatobius j jelskii, T. j. bufo 
(redescribed by LaviLLA, 1988b), T. j. longitarsis and T j walkeri The phenogram (fig. 5) 
of 72 adults including VELLARD's type specimens (MHNSM) shows that: (1) there exist 
seven morphs within the nominal species; (2) these morphs do not coincide with any of the 
subspecies proposed by VELLARD; (3) different morphs inhabit the same locality. Joining 
all specimens to one group causes an error sum of squares of 1.444, that is three times 
greater than in T, culeus or in T. rimac and ten times greater than m T brevirostris. 
Nevertheless, the occurrence of different morphs at all sites — though in different 
frequencies — does not support the validity of the four subspecies proposed by VELLARD. 


Telmatobius rimac Schmidt, 1954 


VELLARD (1955) recognized two subspecies: Telmatobius r. rimac and T. r. meridio- 
nalis. The phenogram (fig. 6) of 42 adults including VELLARD's type specimens shows that: 
(1) there are five morphs; (2) the type material from Ocros and Tupe forms one of two 
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Fig 4 — Phenogram of morphometne similanty between adult Telmatobius culeus which VELLARD 
(1953) assigned to different subspecies. Presentation of data analogous to that in fig. 3. T, c. 
culeus’ Ocama, lave, {sla del Sol, T «. dispar’ Coata, Juhaca; T. c. escomelr Lagunillas; T. c 
exsul Yura; T c fluviatths’ Chucuito, Huayllata, Have; T. c. lacustris. Checayani, Umayo 


main clusters including two mixed morphs; (3) the more recently collected material from 
Obrojillo and Huaytara (situated at the same distance to Ocros in the north and to Tupe 
tn the south) forms the second main cluster and segregates into three more morphs; (4) 
only one individual (Obrojillo 3) directly joins the cluster of the type specimens. This 
unexpected structure of the phenogram does not support the validity of the proposed 
subspecies, but demonstrates suspicious differences between the external morphology of T. 
rimac which were collected forty years ago and those presently inhabiting the streams of 
the Pacific slope of the Andes. 


INTFRSPECIFIC VARIATION 


In order to compare the intraspecific variation with the morphometric differences 
among the nominal species, we performed principal component analysis and discriminant 
analysis on two groups of geographically neighbouring species: (1) the northern Tel- 
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Fig. 5 Phenogram of morphometric similarity between adult Telmatobius jelsku which VELLARD 
(1955) assigned to different subspecies. Presentation of data analogous to that in fig. 3. T. 
jelsku. Acolla; T j bufo Tambo; T. j. longitarsis: Huancavelica, Puguio: T.) walkeri Ayacucho. 
Unclassified T. jelskii were collected in Palian, Cuyrohuas: and Parinacochas 
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Fig. 6. — Phenogram of morphometric similanty between adult Telmatobius rimac which VELLARD 
(1955) assigned to different subspecies. Presentation of data analogous to that in fig. 3 7. r. 


rimac: Ocros, T. r. meridionalis Tupe. Unclassified T rimac were collected in Obroyllo and 
Huaytara 
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matobiinae including Batrachophrynus brachydactylus, B. macrostomus, Telmatobius 
brevirostris and T. carrillae; (2) the southern Telmatobius species T. culeus, T jelskii and 
T. rimac. The main reason for the subdivision of the complete data set on phenotypes into 
groups of three and four species, respectively, was to reduce the number of significant 
canonical variables to two, which permits the distinction of these taxa in two-dimensional 
scatter plots (figs. 7-8). Here, we present only the results based on log,,-transformed data 
because the separation of taxa by discriminant analysis was almost identical in data sets 
of 18 log, -transformed distances and in those consisting of 19 morphometric ratios. 


Generally, the interspecific differences in size (PCI) by far exceeded those 1n shape 
(PC2, PC3). The size effects on PC2 and PC3 appeared to be small, because shearing did 
not notably improve the separation of taxa. Discriminant analysis led to an almost optimal 
separation of species by combining differences in size and shape. 


In the northern group of central Peruvian Telmatobiinae, the first three principal 
components explained 95.0% of the total variance. PCI distinguishes the large B. 
macrostomus from the smailer three species. The plot of PC2 and PC3 scores (fig. 7A) 
shows a wide overlap between B. macrostomus, B. brachydactylus and T. carrillae, whereas 
PC3 distinguishes four of the five T. brevirostris from the other species An almost 
complete separation of the four taxa was obtained by discriminant analysis, only 2 out of 
53 B. brachydactylus were confounded with T. carrillae (fig. 7B. Table III). The separation 
of taxa is mainly based on size (SVL) and interorbital distance (IOD) in CANI, and on 
size (SVL), head shape (HWID, EYE, IOD) and humerus length (HUML) in CAN2. 


In the southern group of Telmatobius species the first three principal components 
accounted for 92.9 % of the total variance. PC1 distinguishes T. culeus from the other two 
species. The plot of principal component scores (fig. 8A) shows a considerable separation 
of T. jelskii from T. rimac based on PC3, but a complete overlap of T. culeus with both 
other species. An almost complete separation of the three taxa was obtained by canonical 
discriminant analysis, only 1 out of 42 T. culeus was confounded with T. jelskii, and | out 
of 72 T. jelskii with T. rimac (fig. 8B, Table IV). The separation of species is mainly based 
on size (SVL, BH) and head shape (ENOSE, ESNOUT) in CANI, and on size (SVL), head 
shape (HWID, ESNOUT) and extremity length (HANDL, TIBL) in CAN2. 


INTERGENERIC VARIATION 


The next step of analysis concerns the morphometric features of each genus and the 
possibility to identify genus-specific morphometric characters or ratios. The data of the 
two Batrachophrynus species form one group, those of the five Telmatobius taxa the other. 
Again, log,,-tranformed data and ratios provided the same degree of group separation A 
highly significant discriminant function was obtained which correctly classifies 80 % of the 
Batrachophrynus and 96 % of the Telmatobius (log, -transformed data, Table V). Distinc- 
tive morphometric ratios (ANOVA, P < 001) were BH/SVL (fig. 9A) and FEML/TIBL 
(fig. 9B). 
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Fig 7. — Plot of (A) principal component scores and (B) discnminant function scores of the four 
northern telmatobiine species: Batrachophrynus brachydactylus, B. macrostomus, Telmatobius 


brevirostris and T. carrillae. Discriminant functions and classification success are given in 
Table HI. 
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southern telmatobuine species Telmatobius jelsku, T. rimac and T. cudeus Discriminant functions 
and classification success are given in Table IV 
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Table Ill. - Discriminant functions to distinguish among Batrachophrynus 
macrostomus, B. brachydactylus, Telmatobius brevirostris and T. carrillae based 
on 18 logio transformed morphometric characters. 


A Statistical significance 


0.985 6135 
0 842 221.5 34 
0.725 5836 16 


B. Unstandardized discriminant function coefficients 


Coefficients 
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C. Classification success 
Predicted group 
Actual group Batrachophrynus | Batrachophrynus | — Telmatobii Telmatobius 
brachydactylus macrostomus brevirostris carrillae 


B, brachydactylus 51 (96%) 0 0 2 (4%) 
B. macrostomus 0 13 (1005) 0 0 

T. brevirostris 0 0 5 (100%) 0 

T. carnillae 0 0 0 53 (100%) 


D Group centroids 


Species 


B. brachydactylus 


B macrostomus 
T. brevirostris 
T. carrillae 
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Table IV. - Discriminant functions to distinguish among Telmatobius culeus, T. jetski 
and T. rimac based on 18 log;o transformed morphometric characters. 


A, Statistical significance 
547 0.919 0.0416 440.3 36 < 0.00001 
271 0.2694 181.7 17 < 0.00001 


Coefficients 
Character 
BH ; 


Constant 


C. Classification success 


Predicted group 


Actual group Telmatobius culeus | Telmatobius jelskit 


Teimatobiųus culeus 41 (98%) 1 (2%) 0 
Telmatobius jelskai 0 71 (99%) 1 (190) 
Telmatobius rimac 0 0 42 (100%) 


D. Group centroids 


Telmatobius culeus 
Telmatobius jelsku 
Telmatobius rimac 
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Table V. Discriminant functions to distinguish the genera Batrachophrynus and 
Telmatobius based on 18 tog), transformed morphometric characters. 


A. Statistical significance 
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C. Classification success 


Predicted group 
Actual group 


Batrachophrynus 53 (80%) 13 (20%) 
Telmatobius 8 (4%) 200 (96%) 


D. Group centroids 


Batrachophrynus -226 
Telmatobius 0.72 
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Fig 9 — Box- and whisker-plot of morphometric ratios which permit the distinction between the 
genera Batrachophrynus and Telmatobius (A) BH/SVL; (B) FEML/TIBL. The dotted line 
represents the means of each genus. B.b.. Batrachophrynus brachydactylus; B.m.: B. macrostomus; w 
T.b.: Telmatobius brevirostris, T ca. T carrillae, Tcu T. culeus, T.j.: T. jelskii, T.r.. T. rimac. e 
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VARIATION BETWEEN ECOTYPES 


Finally, we examined the relationship of morphological traits caused by the 
convergent adaptation to a specific type of habitat (stream/lake) with the phylogenetic 
relationships between the seven taxa. For this purpose, we pooled the data of B. 
macrostomus and T. culeus forming the lake group and compared them with the stream 
group formed by the remaining species. Again, there was no difference between the 
separation of groups based on log,,-transformed data and that based on ratios. The highly 
significant discriminant function correctly classifies 100 % of the stream-inhabitants and 
92 % of the lake-inhabitants (log;,-transformed data, Table VI). Distinctive morphometric 
ratios (ANOVA, P < 0.01) between lake- and stream-inhabitants are HWID/SVL (fig 
10A) and EYE/SVL (fig. 10B). Moreover, SVL of adult lake-inhabitants is considerably 
larger than that of stream-inhabitants (Table I). 


DISCUSSION 


The general similarity 1n síze and shape of the riparian central Peruvian Telmatobii- 
nae does not facilitate a reliable and objective definition of taxa. TRUEB (1979) stated that 
most taxonomic descriptions of Telmatobius are inadequate because many proposed 
diagnostic characters are so subjective that an identification without comparison with type 
material ıs almost impossible. A recent study on the Telmatobius species of northern Perú 
used for the first time multivariate statistics on morphometric data to obtain objective 
criteria for the classification of telmatobiine frogs (Wiens, 1993). We followed this 
approach and successfully applied cluster, principal component and discriminant analyses 
to assess the morphometric variation within and among the central Peruvian taxa of 
Telmatobiinae. Consequently, we offer discriminant functions based on external characters 
which permit the identification of adults with a very low rate of erroneous classifications. 
Moreover, we identified convergent morphological traits distinguishing riparian and 
lake-inhabiting telmatobiine frogs which in the future should be avoided for taxonomic 
conclusions. 


REASSESSMENT OF THE SUBSPECIFIC TAXA 


Within wide ranging species such as T. jelskii and T. rimac we are confronted with the 
problem of interpopulational morphological differentiation, due to the partial geographi- 
cal isolation between different hydrographic systems and valleys. VELLARD (1951, 1953, 
1955) attempted to solve this problem by naming more than half of the populations 
studied up to 1955 at the subspecific level. However, the existence of a complete series of 
intermediate specimens (VELLARD, 1955) between all subspecies emphasizes the more or 
less arbitrary nature of their definition. Nevertheless, due to VELLARD's proposal, the 
segregation of T brevirostris, T. culeus, T jelskii and T. rimac into a total of fifteen (!) 
subspecies is still recognized (Frost, 1985). 


Our attempt to deal with interpopulational variability in taxonomic terms is based on 
the hierarchical grouping of individual phenotypes according to their similarity in external 
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Table VI. - Discriminant function to distinguish between the stream- and lake- 
inhabiting species of the genera Batrachophrynus and Telmatobius based on 18 
logi transformed morphometric characters. 


A. Statistical significance 
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C. Classification success 


Predicted group 
Actual group 


Stream-inhabitant 224 (100%) 
Lake-inhabitant 4 (8%) 


D. Group centroids 


Stream-inhabitant 
Lake-inhabitant 
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Fig 10 — Box- and whisker-plot of morphometric ratios which permit the distinction between 
inhabitants of lakes and of streams (A) HWID/SVL; (B) EYE/SVL. Presentation of data 
analogous to fig. 9. 
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morphology To avoid major biasses due to size, the phenotypes are described by 19 
unweighted morphometric ratios. The phenotypic grouping in a cluster analysis permits an 
objective assessment of similarity between conspecific individuals, even if the sample size 
is small as usually for type material. Principal component and discriminant analyses, in 
contrast, require about 20 cases per predictor which imply sample sizes exceeding by far 
the numbers of type specimens. 


VELLARD’s proposal of subspecies is exclusively based on differences in the external 
morphology of specimens from different localities. Therefore, if the proposed subspecies 
were valid, the phenograms of conspecific individuals collected from different localities (= 
reproductively isolated populations) should show the following structure: (1) specimens 
pertaining to the material originating from one type locality should form a homogeneous 
group (cluster); (2) material assigned to different subspecies should be represented in 
different clusters; (3) conspecific specimens originating from one population without their 
own subspecific status should join as a group one cluster formed by type specimens. 


The phenograms obtained for T. brevirostris, T. culeus, T jelskii and T. rimac do not 
show structures compatible with VELLARD's taxonomic suggestions. The morphometric 
differentiation between the five specimens of T. brevirostris which VELLARD assigned to 
three subspecies is low and the most similar pair of individuals belongs to different 
subspecies (fig. 3). The greatest morphometric difference between the 42 7. culeus exceeds 
three times that between the 7. brevirostris specimens but none of the five morphs 
identified within this data set coincides with any of the proposed subspecies (fig. 4). The 
same applies to the seven and five morphs, respectively, found within the T. jelskii (fig. 5) 
and T. rmac (fig. 6) data sets. As objective morphometric similarity does not correspond 
to VELLARD's subjective grouping of conspecific specimens to subspecies, we conclude that 
the proposed subspecific segregation ts invalid in the four species studied This conclusion 
does not rule out that the definition of subspecies may be useful to describe the speciation 
processes due to the geographic isolation of populations of wide ranging species. However, 
if considered useful, the definition should be based on conceivable, objective criteria better 
than those presented so far. 


What is the meaning of the different morphological traits (morphs = groups of 
conspecific specimens which are morphometrically very similar) which we detected in all 
seven species? As different morphs usually occur at the same locality, i.e. within the same 
population, they probably do not represent taxonomically relevant units. Neither do they 
reflect sexual dimorphism in shape, or size-specific polymorphism. The varying frequencies 
of morphs at different localities indicate that the adaptive value of a morph for the 
increase of individual fitness probably differs among the localities. Thus, intraspecific 
polymorphism may have evolved in response to local environmental factors such as 
altitude, water flow and competition. The occurrence of different morphs at the same 
locality may be due to temporal changes in the local environments. This interpretation is 
supported by the conspicious change of the frequencies of T. jelskii morphs which took 
place in the Mantaro Valley during the last forty years 9 out of 10 specimens collected 
near Acolla in the early fifties (VELLARD, 1955), but only i out of 7 frogs recently collected 
near Palian join the same main cluster, whereas only 1 Acolla specimen joins the remaining 
6 Pahan specimens. Even more impressive 1s the same tendency in T. rimac though 
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referring to geographically distant localities: only one recently collected individual pertains 
to the morphs dominating about forty years before (VELLARD, 1955). Attributing these 
changes to the different duration of preservation seems too simple because older morphs 
still occur (rarely) in contemporary populations, and presently dominating morphs also 
existed (rarely) forty years ago. Instead, there have been dramatic changes in the riparian 
habitats of these species during this period due to the enormous increase in human 
population and the resulting pollution of the rivers and streams. 


In summary, we reject the proposal of subspecies for T. brevirostris, T. culeus, T. 
Jelskii and T. rimac. 


REASSESSMENT OF THE SPECIES STATUS 


The large numbers of studied specimens in most species permit a reliable estimate of 
the intraspecific morphometric variability and the comparison between different species. 
Principal component analysis of the log,,-transformed data did not substantially 
contribute to the resolution between different taxa: the usually large overlap between 
different species reflects their similarity in many aspects of shape. The only exception from 
the rule is T brevirostris which considerably differs from the neighbouring species (fig. 
TA). However, WiENS (1993) obtained a similar low degree of resolution between eight 
Telmatobius species from northern Pert, indicating that PCA is not a powerful tool to 
distinguish between Andean Telmatobiinae. 


In contrast, if groups (— species) are determined a priori, and if the differences 
between them are maximized by canonical discriminant analysis, the same morphometric 
data sets can distinguish between the taxa. It is noteworthy that, despite the bad image of 
morphometric ratios (e.g. BooKsTEIN et al., 1985), discriminant analysis using logo 
transformed data or ratios give very similar results and reach the same classification 
success. The convergent morphological adaptation to similar environmental constraints 
prevents an unequivocal identification at the level of individuals, but the rate of 
erroneously classified individuals is low: 2 B. brachydactylus, | T culeus and 1 T. jelskii 
out of a total of 280 specimens. Wiens (1993) did not provide results of classification 
Success in his analysis, but he also stated that discriminant analysis provided an objective 
base for the distinction between Andean species of Telmatobiinae. In summary, the six 
central Peruvian species as well as the one from southern Perú proved to be well-defined 
taxonomic units which possess external features allowing an objective diagnosis. 


REASSESSMENT OF THE GENERIC STATUS 


LAURENT (1983) recognized three genera within the central Peruvian Telmatobiinae: 
the monotypic genera Batrachophrynus and Lynchophrys, and Telmatobius with four 
species. However, the evidence presented for the change of Batrachophrynus brachydactylus 
to the genus Lynchophrys i$ weak. the main differences from Batrachophrynus macrostomus 
put forward are: (1) shorter third finger; (2) smaller size; (3) the statement that male B. 
macrostomus lack nuptial pads. Points (1) and (2) are convergent morphological 
adaptations to the stream habitat which are shared with all riparian Telmatobius, whereas 
point (3) is simply an error (see figure 16 in SiNscH, 1990). Moreover, both species share 
two morphometric features which are diagnostic for the genus Batrachophrynus despite 
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their adaptation to different habitats: flat body and femur length exceeding tibia length 
Finally, the analysis of allozymes of these species and of three Telmatobius species (SINSCH 
& JURASKE, 1995) clearly demonstrates that B. brachydactylus and B. macrostomus are 
closely related species as originally proposed by PETERS (1873). 


In summary, we do not see any conceivable reason to split the genus Batrachophrynus 
and reject the proposal of Lynchophrys. 


SUMMARY OF TAXONOMIC PROPOSALS 


We summarized in Table VII the reassessment of the taxonomic status of central 
Peruvian Telmatobiinae. Three genera including seven nominal species which segregate 
into 15 subspecies are currently reduced to two genera Batrachophrynus and Telmatobius 
which include seven species without subspecific segregation. 


Finally, we wish to comment on DUELLMAN's (1993) statement that the reading of 
Batrachophrynus macrostomus Peters, 1873 should be changed to B. microstomus. We do 
not agree because there is no doubt about the original naming by PETERS (1873), and in 
all research papers dealing with this species (AviLA Ramon, 1953; CAMARENA, 1953; 
Dunors, 1984, FrELDsA, 1983; GORHAM, 1966; LAURENT, 1983; LaviLLa, 1988a, LYNCH, 
1971, 1978; Macepo, 1950, 1960, 1976, MORALES, 1988; SiNscH, 1986, 1990; VELLARD, 
1951) the original name has been used consequently. 


CONVERGENT MORPHOLOGICAL TRAITS 


Stream habitats require frogs to evolve morphological adaptations which enable them 
to move within the permanent current of water Therefore, it is not surprising that the 
riparian B. brachydactylus shares two diagnostic features with the riparian Telmatobius: 
slim head and large eye diameter. Moreover, adults of all riparian species are smaller-sized 
than those of the lake-inhabiting B. macrostomus and T. culeus. Confusion of these 
convergent morphological adaptations to the same type of habitat with similar morphol- 
ogy due to phylogenetic relationship has led to the creation of the genus Lynchophrys by 
LAURENT (1983). In the Andean Telmatobiinae any taxonomic conclusion based on 
external morphology should be backed up by other kinds of characters because the rate 
of erroneous classification in the morphometric distinction of the genera Batrachophrynus 
and Telmatobius 1s considerably greater than that in the distinction between stream- and 
lake-inhabitants: 20 % versus 8 %. Thus, the effect of convergent lines of development 1s 
probably great in this group of frogs. 


GEOGRAPHICAL DISTRIBUTION 


The geographical range (fig. 2) of the six central Peruvian species is still relatively 
unknown. We know little about the northern extension of the ranges of T. brevirostris, T. 
carrillae and T. rimac, and about the southern range of T. jelskii. The northern gap 
between the species surveyed by Wrens (1993) and those in this study 1s subject to a recent 
study (SALAS, in prep.). Further attention should be paid to the exact limits of distribution 
of the riparian B. brachydactylus in relation to those of the neighbouring Telmatobius 
species. We do not even know if B. brachydactylus and T. jelskit, both present in the 
hydrographic system of the Mantaro river, can coexist at the same locahty. 
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Table VII, - Alphabetical synonymy of the telmatobiine species revised 1n this study. 


Lynchophrys brachydactyla (Peters, 1873) 
Batrachophrynus macrostomus Peters, 1873 
Telmatobius brevirostris Vellard, 1955 


Batrachophrynus brachydactylus Peters, 1873 
Batrachophrynus macrostomus Peters, 1873 


revirostris Vellard, 1955 
arvulus Vellard, 1955 
unctatus Vellard, 1955 


Telmatobius brevirostris Vellard, 1955 
Telmatobius brevirostris Vellard, 1955 
Telmatobius brevirostris Vellard, 1955 


Telmatobius carrillae Morales, 1988 Telmatobius carrillae Morales, 1988 


Teimatobius culeus (Garman, 1875) 
T. culeus cuteus (Garman, 1875) 
T. culeus dispar Vellard, 1953 
T. culeus escomeli Angel, 1923 
T. culeus exsul Vellard, 1951 
T culeus fluviatilis Vellard, 1953 
T culeus lacustris Vellard, 1953 


Telmatobius culeus (Garman, 1875) 
Telmatobius culeus (Garman, 1875) 
Telmatobius culeus (Garman, 1875) 
Teimatobius culeus (Garman, 1875) 
Telmatobius culeus (Garman, 1875) 
Telmatobius culeus (Garman, 1875) 
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Telmatobius jelskii (Peters, 1873) 


Telmatobius rumac Schmidt, 1954 


T. jelskit jelski (Peters, 1873) 

T. jelskii bufo Vellard, 1955 

T. jelskii longitarsis Vellard, 1955 
T. jelskii walkeri (Shreve, 1941) 


T. rimac rimac Schmidt, 1954 
T rimac meridionalis Vellard, 1955 


Telmatobius jelskii (Peters, 1873) 
Teimatobius jelskii (Peters, 1873) 
Telmatobius jelskii (Peters, 1873) 
Teimatobius jelskii (Peters, 1873) 


Telmatobius rimac Schmidt, 1954 
Telmatobius rimac Schmidt, 1954 
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Nevertheless, there can be little doubt that the genus Batrachophrynus is endemic to 
central Perú, as already mentioned in the original description and later in VELLARD (1951). 
The comment on distribution in Frost (1985) — “Andes of southern Peru and Bolivia" 
— is obviously wrong, as is the citation of “Lynchophrys brachydactyla" occurring in the 
northern Andes of Bolivia (DE La Riva, 1990). 


CONCLUSION 


At this stage of morphometric analysis we refrain from phylogenetic considerations 
because of the convergent adaptations of the external morphology to the same type of 
habitat. In the next step of analysis, we use allozyme variation within and between 
telmatobiine species to approach the phylogenetic relationships and to compare them with 
phenetic relationships based on morphometry (SinscH & JuRASKE, 1995). 


In conclusion, this paper provides an objective, diagnostic method to assign central 
Peruvian Telmatobiinae to the presently known species, based exclusively on external 
characters which are easy to measure. Thus, multivariate statistics, specifically canonical 
discriminant analysis, have proven again to be an useful tool in the classification of 
amphibians. 


RESUMEN 


Se revisa el estado taxonómico de los Telmatobiinae de la región central del Perú, en 
base a la variación intraespecífica de 18 parametros morfométricos que presentan las seis 
especies (con nueve subespecies) reconocidas actualmente, las cuales se agrupan en tres 
géneros (Batrachophrynus, Lynchophrys y Telmatobius). Aphcando los análises de 
componentes principales y de discriminación se reconocen dos géneros (Batrachophrynus 
y Telmatobws) incluyendo seis especies (Batrachophrynus brachydactylus, B. macrostomus, 
Telmatobius brevirostris, T. carrillae, T. jelski, T. rmac) y mnguna subespecie. Dos 
caracteres diagnósticos y externos distinguen las especies de Batrachophrynus y de 
Telmatobius, otros dos son considerados adaptaciones convergentes que diferencian a los 
habitantes de arroyos de aquellos que habitan las lagunas 
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APPENDIX I 
SPECIMENS EXAMINED 
Batrachophrynus brachydactylus 


Perú: (1) Departamento Cerro de Pasco: Caza Pato, 9 males, 3 females, MHNSM 
1311, 1314, 1316, 1319, 1324-1325, 1331-1336; (2) Departamento Junin: brook near 
Ondores: (a) first sample, 9 males, 10 females, URP 054-072; (b) second sample, 11 males, 
11 females, collected by A. SALAs, March 26, 1992. 


Batrachophrynus macrostomus 
Perú: Departamento Junin: Junín Lake, Carhuamayo, 6 males, 7 females, 8 juveniles, 
collected by M. ANTIGNANI, February 20, 1992. 
Telmatobius brevirostris 


Perú: Departamento Huanuco: (1) Ambo, Chasqui, 1 male, 1 female, MHNSM 3736, 
7676 (syntypes of T. b. brevirostris), (2) Ambo, Caina, 2 females, MHNSM 7666-7667 
(syntypes of T. b. parvulus), (3) Santa Maria del Valle, 1 male, MHNSM 7681 (holotype 
of T. b. punctatus). 


Teimatobius carrillae 


Perv: Departamento Ancash: (1) Yuracyacu: (a) first sample, 4 males, 2 females, 
MHNSM 1528 (holotype), 1544-1545, 3932-3934 (paratypes); (b) second sample, 9 males, 
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12 females, URP 001-021; (2) Huikia: (a) first sample, 5 males, 1 female, MHNSM 
6681-6687; (b) second sample, 6 males, 2 females, URP 022-029; (3) Huaychopampa, 8 
males, 4 females, URP 030-041. 


Telmatobius culeus 


BoutviA: Lake Titicaca, Isla del Sol, 1 male, 1 female, MHNSM 7769-7770 (assigned 
to T. c. culeus). 


PERU: (1) Departamento Arequipa: Yura, Arequipa, 1 female, MHNSM 7678 
(syntype of T. c. exsul); (2) Departamento Puno: (a) Azangaro, Hacienda Checayani, | 
male, | female, MHNSM 7673-7674 (syntypes of T. c. lacustris), (b) Lake Lagunillas, 6 
males, 4 females, MHNSM 7768, 7776-7777, 7785-7786, 7806-7807, 7823-7825 (assigned to 
T. c. escomeli); (c) Lake Titicaca, Ocama, 3 males, 3 females, MHNSM 7779-7784 
(assigned to T. c. culeus); (d) Rio Coata, Juliaca, 2 males, MHNSM 7771-7772 (syntypes 
of T. c. dispar); (e) Rio Ilave, Chucuito, 1 male, 1 female, MHNSM 7754-7755 (syntypes 
of T. c. fluviatilis); (f) Rio Ilave, Huayllata, 1 male, 2 females, MHNSM 7812-7814; (g) 
mouth of Rio Ilave, 2 females, MHNSM 7773-7774 (assigned to T. c. culeus), and 3 
females, MHNSM 7812-7814 (assigned to T. c. fluviatilis); (h) Rio Juliaca, Puno, 4 males, 
4 females, MHNSM 7766-7767, 7787-7789, 7793-7794 (assigned to T. c. dispar); (i) 
Umayo, 34 km NO of Puno, | female, MHNSM 7811 (assigned to T. c. lacustris). 


Telmatobius jelskii 


Perú: (1) Departamento Ayacucho: (a) Ayacucho, 4 males, 6 females, MHNSM 
12202, 12206, 12213-12214, 12217, 12219-12220, 12222, 12225, 12899 (assigned to 7T. j. 
walkeri), (b) Parinacochas, 6 males, 6 females, MHNSM 12838, 12841, 12883, 12901, 
12904-12910; (c) Puquio, 2 males, | female, MHNSM 7642-7643, 7645 (syntypes of T. j. 
longitarsis); (d) Tambo, 2 males, 4 females, MHNSM 7646-7651 (syntypes of T. jelskii 
bufo), (2) Departamento Huancavelica: Huancavelica, 4 males, 1 female, 2 juveniles, 
MHNSM 7639-7641, 7660-7661, 7663-7664; (syntypes of T. j. longitarsis); (3) Departa- 
mento Junin: (a) Huancayo, Acolla, 3 males, 7 females, MHNSM 6903-6906, 6909-6914 
(assigned to T. j. jelskii); (b) Huancayo, Palian, Rio Shullcas, 5 males, 2 females, 14 
juveniles, collected by U. SiNscH & V. CANALES, February 19-25, 1992; (c) Tarma, 
Cuyrohuasi: (i) first sample, 5 males, 5 females, URP 90066-90075; (ii) second sample, 6 
males, 3 females, collected by A. SALAS, March 24, 1992. 


Telmatobius rimac 


Perv: (1) Departamento Ancash: Ocros, 7 males, 3 females, MHNSM 6935-6936, 
6941-6942, 6944-6945, 6950-6951, 6953-6954 (assigned to T. r. rimac); (2) Departamento 
Lima: (a) Canta, Obrojillo, Rio Chillon, 8 males, 7 females, 7 juveniles, collected by J. 
IcocHEA, March 8-9, 1992; (b) Canta, Quebrada Huaytara, 9 males, 4 females, 2 juveniles, 
collected by J. ICOCHEA, March 8-9, 1992; (c) Tupe, 2 males, 2 females, MHNSM 
7656-7659 (syntypes of T. r. meridionalis). 
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